Scattering of polarized laser light by an atomic gas in free space: a QSDE approach 
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We propose a model, based on a quantum stochastic differential equation (QSDE), to describe 
the scattering of polarized laser light by an atomic gas. The gauge terms in the QSDE account for 
the direct scattering of the laser light into different field channels. Once the model has been set, we 
can rigorously derive quantum filtering equations for balanced polarimetry and homodyne detection 
experiments, study the statistics of output processes and investigate a strong driving, weak coupling 
limit. 
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I. INTRODUCTION 

Many recent experimental 0, 0j S 0] and theoretical 
@j B 0, B [I, 03 works have been based on a simple ex- 
perimental scenario, in which a polarized atomic gas is 
continuously probed with a polarized off-resonant opti- 
cal beam (Fig. [T]). By measuring the Faraday rotation 
of the optical polarization resulting from the interaction, 
one can in principle prepare conditionally spin-squeezed 
states or perform quantum metrology tasks, e.g. estimat- 
ing a magnetic field that rotates the spins. 

Central to the description of these experiments is the 
quantum filtering equation, which propagates the expec- 
tation value of the atomic gas observables conditioned on 
prior measurement results. The conditional expectation 
is the mean least squares estimate of an atomic gas ob- 
servable given the observations thus far. The conditional 
expectations of 'all' atomic observables can be summa- 
rized in an information state 7r t . The filtering equation 
propagates this information state in real time. 

In quantum optics the filtering equation is often re- 
ferred to as the stochastic master equation (llj . For the 
polarimetry example considered here, previous modelling 
efforts have either produced an unconditional description 
Q or arrived at a conditional description by heuristically 
'adding the usual measurement terms' [i[ in analogy with 
a physicallly different homodyne measurement scheme 
with only a single polarization mode @, 0| ■ I n this arti- 
cle we treat the conditional evolution of the state (due 
to detection of Faraday rotation with a polarimeter) in 
a rigorous manner, allowing the atomic system to medi- 
ate exchange between two orthogonal optical polarization 
modes. In particular, we derive the quantum filtering 
equation from an underlying quantum stochastic model, 
i.e. the quantum stochastic differential equation (QSDE) 
governing the interaction of the atomic gas with the laser 
light. 

Formal quantum filtering theory was pioneered by V.P. 
Belavkin in [l2|, [U[ using martingale techniques (see also 
[IH). We here employ the reference probability method, 
based on the quantum Bayes formula fl5l . Hfl. to obtain 
the quantum filter from the QSDE (see also |17j)- 

The QSDE model we use here is based on a simple 
Faraday Hamiltonian, H = nF z S z , where k is a small 



interaction strength prefactor, S z is a Stokes operator 
measuring the circularity of optical polarization and F z 
is the z-componcnt of the collective atomic spin. Under 
this Hamiltonian, photons with a right circular polariza- 
tion rotate the collective atomic spin over a positive angle 
« along the z-axis, while photons with a left circular po- 
larization rotate the collective spin over a negative angle 
— ft. With linearly polarized light, the angle of linear po- 
larization will Faraday rotate by a degree proportional to 
the z-component of the spin. Note that we entirely ne- 
glect 'tensor' terms of the interaction Hamiltonian (non- 
linear in individual spin operators) which arc important 
near resonance with realistic atoms of spin greater than 
1/2 [|. We have also omitted the evolution due to any 
driving magnetic field, e.g. H = jBF y , purely for rea- 
sons of simplicity, it can easily be added at the end. 

In our QSDE-description, the Faraday interaction is 
described as a 'direct' scattering process, without coher- 
ent absorption and re-emission. This is a consequence of 
the fact that the interaction Hamiltonian is derived from 
an approximation in which the excited states are adiabat- 
ically eliminated At present, however, no mathemat- 
ically rigorous treatment of this elimination is available 
in the literature (see [l8| for rigorous results on the adi- 
abatic elimination of a leaky cavity mode). Therefore, 
we have chosen to directly base our QSDE model on the 
Faraday Hamiltonian without proceeding through a rig- 
orous Markov limit 0, [2(| followed by adiabatic elimi- 
nation of the excited states. Mathematically, the direct 
scattering is represented by gauge-terms in the QSDE 

Having set the underlying model, i.e. the QSDE, we 
rigorously derive the quantum filtering equation for the 
balanced polarimetry setup and for homodyne detection 
of the y-polarized channel. We investigate the statistics 
of the output processes for these two experiments and 
take a limit where the driving laser power a 2 goes to in- 
finity but where the product M = K 2 a 2 is kept constant 
(k is the parameter that couples the field to the atomic 
gas). We show that in this strong driving, weak cou- 
pling limit the statistics of the output processes for the 
balanced polarimetry experiment and the homodyne de- 
tection experiment are equivalent. Furthermore, we show 
that in the strong driving, weak coupling limit we obtain 
the quantum filter that has already been intuitively as- 
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sumcd in the literature [4|]. 

The remainder of this article is organized as follows. 
Section [TT] introduces the fundamental noises and the 
quantum stochastic calculus, and section IIIII sets our 
QSDE model. Section HVl derives the filter when counting 
in the 45 degrees rotated sy-basis (balanced polarimc- 
try), and Section [V] derives the quantum filter for the 
homodyne detection experiment. In sections I VII and IVIII 
we study the statistics of the observation processes, and 
investigate the strong driving, weak coupling limit. We 
close the paper with a discussion of the results obtained. 



II. THE QUANTUM CALCULUS 

One polarized photon in a beam of light can be de- 
scribed by the one particle space 

H = C 2 (g>L 2 (R) L 2 (R;C 2 ), 

of C 2 -valued quadratically integrable functions on the 
real line. The polarized light field is described by the 
bosonic Fock space J~(7~l) over Ti. 



A\ and A l t * are called the annihilation and creation pro- 
cesses, respectively. The processes Af are called gauge 
processes. Formally, we can write the noises as A\ = 
/ a l s ds, A" = J Q a"ds and Af = J Q a"a{ds where a\ 
and a? s are the usual Bose fields. Mathematically, the 
objects a\ and a{ are ill-defined and therefore we resort 
to the definition of Eq. ([1]). The formal expressions do 
show very explicitly though, that the operator Af counts 
the number of photons with a polarization in the e% di- 
rection up to time t and that the operator Af scatters 
the polarization of a photon from the ej direction to the 
direction. 

We will usually work in the basis {e x , e y } which phys- 
ically corresponds to an orthonormal basis in the plane 
orthogonal to the direction of propagation of the light. 
Apart from this basis, we also use the circular basis 



given by {e + 



le 



;)/>/2), 



and the 45 degrees rotated xy-basis given by {e^ = 
(e x + e y )/\/2, e n = (e x — e y )/V2}. Given the defini- 
tions in Eq. |T]) it is easy to work out how the noises 
transform under basis transformations. For example, we 
have 
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n=l 



which enables arbitrary superpositions between states 
with a different number of photons. Note that photons 
are bosons and therefore need to be described by sym- 
metric wavefunctions. For an / £ Ti we can define the 
exponential vector e(f) in T(Ti) by 



e(/) 
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We call the span of the exponential vectors the exponen- 
tial domain. The exponential domain is a dense set in 
T(TL) and we allow ourselves the freedom to only pro- 
vide the definition of the fundamental noises (a little 
further below) on this domain. If we normalize the ex- 
ponential vectors then we obtain the coherent vectors 
ip(f) = cxp(— i||/|| 2 )e(/). An important vector is the 
vacuum vector, given by $ = tp(0) = e(0) = 1 © © . . .. 
The vacuum state <f> = • $) is obtained by taking 
inner products with the vacuum vector. 

If we choose an orthonormal basis {ei, e-i\ in C then 
we can decompose every / € L 2 (R;C 2 ) along this basis, 
i.e. / = f%ei + $2^2 with f% and /2 in L 2 (R). We now 
introduce the fundamental noises _Al, A\* and Af on the 
exponential domain by (see also [2~U, I22I, |23| ) 



mi) = / fi(s)ds e(f), 



(e(g),Afe(f)) = Ug^ds) (e(g),e(f)), (1) 
(e(g),A?e(f))= ^gMfA^ (e(j),e(/)). 



A^ + = ^(Ar 



Af 



0- 



Denote by f) the Hilbert space of the atomic gas. The 
space of the combined system of atomic gas and field 
together is then given by Fj <g> !F(H). Define H t ] — C 2 © 
L 2 (-oo, t] and H[ t — C 2 © L 2 [t, 00). For all t the bosonic 
Fock space splits in a natural way as a tensor product 
F{U) = F{H t \) © T{H[ t ). A process L s , (s > 0) on 
I) © T^H) is called adapted if L s acts nontrivially only on 
t) © !F(H S ]) and is the identity on T{TL^ S ) for all s > 0. 

Hudson and Parthasarathy defined stochastic in- 
tegrals of adapted processes L s against the fundamen- 
tal noises, i.e. they gave meaning to the expression 
X t = X + j Q L s dM s where M s is one of the fundamental 
noises A\ , ^4" or A l J . The expression can be written in 
shorthand as dX t = L t dM t . More importantly, Hudson 
and Parthasarathy [13] provided the calculus with which 
these stochastic integrals can be manipulated in calcula- 
tions. The calculus consists of the following. Suppose X t 
and Yt are stochastic integrals, i.e. dX t = L\dM\ and 
dY t = L\dM 2 where L 1 and L 2 are adapted processes 
and M 1 and M 2 are fundamental noises, then the prod- 
uct X t Y t is itself a stochastic integral. Moreover, the 
product X t Y t satisfies the following quantum Ito rule, 
(partial integration rule) 

d(X t Y t ) = X t dY t + (dX t )Y t + dX t dY u 

where to evaluate dX t dY t we use that the increment dM t 
of a fundamental noise commutes with all adapted pro- 
cesses, and products dM}dM 2 are given by the following 
quantum Ito table 22 1 
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dM l \dM 2 


dAi* dk\ j dA\ 


dAf 





dkf 


MA** <MA? 


dA\ 


S ki dt 5 ki dAl 



and all products dM t dt and dtdM t are zero. As an ex- 
ample, suppose dX t = L\dA\ and dY t = L 2 dA\* , then 
d(X t Y t ) = X t L\dA? + LlYdAl + L\L\dt. 

It can be shown [2(| that in the weak coupling 
limit (a Markov limit) QED models converge to quantum 
stochastic models, i.e. in the limit the unitary time evo- 
lution U t satisfies a quantum stochastic differential equa- 
tion in the sense of Hudson and Parthasarathy. Usually, 
the QSDE obtained via the weak coupling limit can be 
simplified further by adiabatic elimination of degrees of 
freedom of the initial system. See for instance [l8| for rig- 
orous results on the adiabatic elimination of a leaky cav- 
ity. We, however, are interested in the adiabatic elimina- 
tion of the excited states of the atoms in the atomic gas. 
Unfortunately, at present, no rigourous results on this 
kind of adiabatic elimination are available. Therefore we 
choose not to go through a weak coupling limit / adiabatic 
elimination procedure here, but rather write down a phc- 
nomenological QSDE based on the Faraday interaction, 
see Eqs. © and © below. See Q for a derivation of the 
Faraday Hamiltonian of Eq. (|2|) via usual non-rigorous 
adiabatic elimination methods. 



III. THE MODEL 

The interaction between the laser light and the spin 
polarized atomic gas is governed by the Faraday interac- 
tion given by 

Hdt = 2nF z S z dt = nF z (dA+ + - dAr~ ) . (2) 

Here k is a coupling parameter, F z is the z-component 
of the collective spin vector of the atoms, and 2S Z = 
a^*af — aF*ar is the z-component of the stokes vector S 
of the polarized light. The time evolution of the coupled 
system of light and atomic gas together is given by the 
exponential (the superscript distinguishes C/ t ° from Ut 
to be introduced later) 

U° = cxp (ij nF z (dk++ -dA~-)j. 

Since A^~ + and are jump processes, the Ito rule leads 
to a quantum stochastic differential equation (QSDE) 
[22I ] that contains the following difference terms (U§ = I) 

dU? = {(e"^ - l)dA++ + (« " ;/ - l)dA t — } U°. 

(3) 

That is, right circular polarized photons rotate the col- 
lective spin of the atoms over an angle k along the z-axis, 
whereas left circular polarized photons rotate the collec- 
tive spin of the atoms over an angle — k along the z-axis. 




Probe 
laser 



FIG. 1: Schematic depicting balanced polarimetric detection 
of laser light after interacting with a polarized cloud of atomic 
spins via the Faraday Hamiltonian. The light is initially lin- 
early polarized along the x direction. After the interaction, 
the light carries off information about the atomic gas encoded 
in a small optical polarization rotation. The light is measured 
in the £-n basis rotated 45 degrees from the x-y basis, such 
that without the atomic gas the mean output of the polarime- 
ter is balanced to zero. The change of measurement basis is 
achieved with the waveplate located just before the polarizing 
beamsplitter. 



If we express the gauge processes in the linearly polarized 
xy-basis, then Eq. ^ reads (11® = I) 

dU° = { ( cos(kF 2 ) - 1) (dAf + dXf y ) - 

sm( K F z )(dA* v -dAf)}<7 t . 

(4) 

The second term shows that the interaction can scatter re- 
polarized photons to y-polarized photons and vice versa. 

Initially, the atomic gas is in an a;-spin polarized state, 
denoted p, and the field is in an x-polarized coherent state 
ip x {f) which represents the driving laser. The function 
/ G L 2 (M. + ) gives the phase and amplitude of the driving 
laser field at every time t G K + . In computations it 
is often convenient to work with respect to the vacuum 
state 4> = ($, • $) for the field. We can obtain a coherent 
state by acting with a displacement or Weyl operator 
W x (f) on the vacuum vector 

If we work with respect to the vacuum state r/>, then we 
have to sandwich all operators with the Weyl operator 
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W x (f). An observable S of the combined system of 
atomic gas and field up to time t is therefore at time 
t given by 

MS) = w x (fyutsu?w x (f) 
= w x (f t yu° t *su°w x (f t ). 

Here, ft denotes the function / truncated at time f, i.e. 
f t (s) = f(s) for all s < t and f t (s) = for all s > t. 
The relation Eq. ([5|) follows since all the operators split 
as a tensor product at time t and U® and S act as the 
identity operator after time t. Since W x {f) is unitary, it 
then cancels against its adjoint for the part that is after 
time t. 

It can be shown [23| that W x (ft) satisfies the following 
QSDE (W x (f a )=I) 

dW x (f t ) = {f{t)dA x t * ~l(t)dA x - \\f{t)\dt}w x {f t ). 

Defining Ut = U®W X (ft) and using the quantum Ito rule 
H2], we obtain (U = I) 

dU t = {(cos(kF z ) - l)(dA xx + dAf) - 
sm( K F z )(dA xy - dAf) + /(f) cos(nF z )dA x * - (6) 
J(t)dA x + f(t)sm(nF z )dA?* - ±\f(t)\ 2 dt}u t . 

Summarizing, we work in the state P := p ® (j>, the time 
evolution of (adapted) observables S is given by jt(S) = 
U*SU t , with U t given by Eq. ©. 

IV. THE QUANTUM FILTER 

After the interaction, the light carries off information 
about the atomic gas. Therefore, measuring the field 
will enable us to make inference about the atomic gas 
observables. Let us suppose that we are counting the 
photons with a polarization along the = l/\Z2(e x + e y ) 
axis, and that we are separately counting the photons 
with a polarization along the e v = l/y/2(e x — e y ) axis, 
see Fig. [T] That is, our observations are given by 

Yi = u;k?u t = \u;{a xx + Af + Af + K\ x )U t , 

1 (7) 

Y? - UfAfUt = -U;(k x t x + Af - A xy - Af)U t . 

Let X be an atomic gas operator, its time evolution is 
given by 

j t (x) - u;xu t . (8) 

Eq. |jgj is called the system. Together Eqs. © and Q 
form a system-observations pair. 

It is easily checked that [Y t a 7 Yf] = for all a,/3 G 
{^, 77} and for all t, s > 0. This is called the self- 
nondemolition property and ensures that our observa- 
tions are simultaneously observable classical processes. 



Furthermore, it can be shown that [jt(X),Y"} = for 
alH > s > and a € {^, rj}. This is called the nondemo- 
lition property. Together the self-nondcmolition and the 
nondemolition property ensure the existence of the con- 
ditional expectation ¥(j t (X)\y t ) of a system operator at 
time t on the observations up to time t. Since the con- 
ditional expectation is linear in the atomic gas operators 
X, we can define an information state ir t on the atomic 
gas system by 

n t (X)=¥(jt{X)\y t ). 

Note that irt is a stochastic state since it depends on the 
observations and Y v up to time t. 

It is the goal of quantum filtering theory to obtain a 
recursive stochastic differential equation that propagates 
the information state 7r t in time. Our approach here is 
based on the reference probability method [HI, In 
these references, the interested reader can find further 
details on the exposition below. 

Our first step is one of mere convenience. It is a change 
of picture that will simplify subsequent calculations. Let 
Wt be given by Wo = I and 

dW t = \j(t)dA x - f(t)dA x * - i|/(*)|dfc| W t . 

Note that W t is the adjoint of W x (f t ). Now define U' t = 
W t U t , where U t is given by Eq. ©. It easily follows from 
the quantum Ito rule [22j that 

dU't = { ( cos(kF z ) - 1) (dA xx + dAf) - 
sin(KF z )(dA xy - dA yx ) + f(t)(cos{ K F z ) - l)dA x * - 
J(t)(cos{nF z ) - l)dA x + f(t)sm(nF z )dA y t * - 

/(f) sm(KF z )dA y t + |/(f)| 2 (cos( K F z ) - l)df}f/ t '. 

(9) 

Define a new state on the combined system of atomic gas 
and field by Q^S) = F(UfSUj.). To complete our change 
of picture we need to sandwich the observables with the 
opposite rotation. That means that the system Eq. ([8]) 
is now simply given by U' t U;XUtU' t * = W t XW* = X. 
In the last step we used that X acts on the atoms and 
is the identity on the field and vice versa for Wt- In the 
new picture the observations read 

zl = upfui* = w t Af w;, 

= U' t Y?U' t * = WtAfW*. 
Using the quantum Ito rule, it easily follows that 

dZ\ = dAf + \{7(t){dA x + dA y ) + 

f(t)(dA x *+dAr) + \f(t)\ 2 dt), 

dZ? = dAT + \(j{t){dA*-dA y ) + 

f(t)(dA x *-dAD + \f(t)\ 2 dt). 



5 



Denote Cf = U{.ytUt* , i.e. Ct consists of the processes 
and Z v up to time t. It can easily be shown that the con- 
ditional expectations in the two different pictures are re- 
lated by F(j t (X)\y t ) = Ul*Q*(X\C t )Ul. That completes 
our discussion of the change of picture. We will now focus 
on deriving an equation that propagates Q t (X\Ct). 

At the heart of the reference probability method is the 
following quantum Bayes formula flU . Il6j . Let V be an 
operator that commutes with Zf for all a G an d 
for all < s < t. Moreover, suppose that V*V > 
and that ¥(V*V) = 1. Then we can define a state Q by 
Q(S) = P(V*SV), and for all operators X that commute 
with Zf (a G {£, ?/}, < s < t), we have (see [H [H] for 
a proof) 

_ P(U*XU|C t ) 



P(U*U|C t ) ' 

We would like to apply the quantum Bayes formula to 
Q*, i.e. with V = U' t . However, Eq. © shows that U' t 
is driven by noises that do not commute with Z" (a G 
{Sj^IjO < s < t), i.e. U[ itself does not commute with 
the Z^s. 

The following trick [13] solves this problem. Suppose 
V{ satisfies the QSDE 

dVt = {(cos(kF z ) + sin{KF z ) - l)dZ^ + 

(cos(kF z ) - sm( K F z ) - l)dZ?\vJ. 

(11) 

Then, the coefficients of dA**, dA\* and dt are the same 
as in Eq. ©. Since dAf and dA"^ (a, (3 G {x,y}) are 
zero when acting on the vacuum vector we therefore 
have that for all operators S [l7| 

Q^S) = F(ul*sul) = F(yfSVj). 

Moreover, since V{ is driven by z\ and Z%, it com- 
mutes with C t , and we can therefore apply the Bayes 
formula with V = V( . That is, summarizing what we 
have achieved thus far 

(12) 

The next step is to find the equation that propagates 
P(V t '* XV t '\Ct) in time. Using the quantum Ito rule we 
find 



dV{*XV{ = Vf{LtXlt - X)V t 'dZt 



with 



vj*(pxif> - X)V{dZ r >, 



L> = cos(nF z ) + sin(/s-F*), 
L v = cos(kF z ) — sin(/c.F z ). 



(13) 



(14) 



We can write Eq. (|13p in integral form and approximate 
the stochastic integrals in the usual way with simple pro- 
cesses. If we proceed by taking the conditional expec- 
tation P( • \Ct), then we can pull the integrators which 
are elements of Ct out of the expectation. Furthermore, 
the conditional expectation ¥(L s \Ct), (0 < s < t) of an 
adapted process L equals ¥(L S \C S ). In this way we obtain 

dF{Vl* XVl\C t ) = V(Vl*(L^XL" - X)V(\C t )dZl + 
¥{Vl*{L^XL^ - X)Vl\Ct)dZ%. 

Now define a t (X) = U'*¥{V(* XVl\C t )U' t for all atomic 
operators X . Using the quantum Ito rule, we obtain the 
linear version of the quantum filtering equation 



do- t {X)=a t (C(X))dt 



a t (L^XL^ - X) ( dY$ 
aAL^XU 1 - X 



dY? 



\f(t)\ 2 dt 
\f(t)\ 2 dt 



(15) 



where the Lindblad generator C is given by 
C(X) = \ f(t)\ 2 (sm(KF z )Xsm{ K F z ) + 

cos(kF z )X cos(kF z ) - x) , 



(16) 



for all atomic operators X. Now recall from Eq. (1121) that 
TTt{X) = o~t(X)/at(I), which is a quantum version of the 
classical Kallianpur-Striebel formula. Using the Ito rule 
once more, we obtain the following quantum filter 

d-K t (X)=TT t (C(X))dt + 

^§w-" ix ) (^->i 2 *.<^)+ 

(^T"H (^->)IM^)4 

The processes dY? - ^\f(t)\ 2 n t (L^L^)dt and dY? - 
^\f(t)\ 2 Tr t {L ri L T, ^dt are called the innovations or inno- 
vating martingales. It can indeed be shown [l3l . [lij that 
the innovations are martingales with respect to the fil- 
tration y*t and the measure induced by P. 



V. A DIFFERENT SETUP: HOMODYNE 
DETECTION 

For the remainder of the paper we assume that f(t) = 
ae"*' with a real and (j>t in [0,27r). Now suppose that 
instead of the balanced polarimetry setup described in 
the previous section, we use a homodyne detection setup 
to measure the ^-component of the output light, sec Fig. 
[2] It is well known [ill, 0] that for such a homodyne 
detection setup the observations are given by 



Y t = Ulie-^AX + e***A%*)U t . 



(17) 



6 



That is, for homodyne detection of the y-channel the 
system-observations pair is given by Eqs. ||SJ) and (fTT|) . It 
is easily checked that the homodyne system-observations 
pair satisfies the self-nondemolition and nondemolition 
properties, meaning that the conditional expectation 
P(jt(X)\yt) is well-defined. Here 3^t denotes the homo- 
dyne observations of Eq. (|17[) from time up to time 
t. Wc will now derive the filter for the corresponding 
information state n t (X) = ¥(j t {X)\y t ). 

Our first step is again one of convenience. We change 
to the Schrodinger picture by defining the following state 
on the combined system of atomic gas and field together 
Q*(5) = P(U;SU t ). In the Schrodinger picture our sys- 
tem is simply given by Utjt{X)U^ = X and the observa- 
tions are given by 



U t Y t U t * 



e^Af 



Denote Ct = [/t^tl^*, i.e. Ct consists of the process Z up 
to time t. It can easily be shown that the conditional 
expectations in the Heisenberg and Schrodinger pictures 
arc related by P(j t (X)\y t ) = U^{X\C t )U t . 

To compute Q t (X\Ct) wc would like to use the Baycs 
formula. Suppose V t satisfies the following QSDE (V = 

I) 

2 

2 (18) 

Then, the coefficients of dA x *, dA\* and dt are the same 
as in Eq. ©. Therefore we have Q\S) = P(U?SU t ) = 
P{V t *SV t ). The equation for V t is driven by Z and A x * , 
both commute with C t , i.e. Vt commutes with Ct- That 
means we can now apply Bayes formula with V = Vt to 
obtain 

p( Jt (x)\y t ) = u t *q\x\c t )u t = 

\.jty >\st) S£ v Ut*P(V t *Vt\Ct)Ut 

Using the quantum Ito rule we find 
dVtXVt = V t *C(X)V t dt + 

ae-^V* cos{nF z )XV t dA x t + 
ae^VfX cos(KF z )V t dAf + 
aV t * (sm(nF z )X + X sm(nF z ))VtdZt. 

where C is given by Eq. (TTB]). Since dA x and dAf* 
are independent of Ct and since vacuum expectations of 
stochastic integrals with respect to dAf and dA x * are 
zero, we find in an analogous way as before 

dP(V t *XVt\C t ) =P(y t *C{X)V t \Ct)dt + 

aP(y t *(sin{KF z )X + X sia( K F z ))V t C^jdZ t . 

Now introduce a t (X) = U;P{V* XV t \C t )U t for all atomic 
gas operators X . Using the quantum Ito rule, wc obtain 
the linear homodyne filtering equation 



da t (X) =a t (C{X))dt 




Ignore x-component 



50/50 



Local Oscillator 
(y-polarized) 



Polarizing 
Beamsplitter 



Post- Interaction 
Probe Beam 



FIG. 2: In the homodyne detection setup, the detection ap- 
paratus in the dashed box of Fig. [T] should be replaced with 
the apparatus depicted schematically here. As in Fig. [1] the 
light is initially polarized along x and the polarization rotates 
slightly due to the interaction with the atoms. Here, however, 
the strong x component is split off and ignored while the weak 
y component is sent to a standard homodyne setup. The y 
component is mixed at a 50/50 (non-polarizing) beamsplit- 
ter with a strong local oscillator beam also polarized along 
the y direction and derived from the same laser as the probe 
beam. The photocurrent representing the interference signal 
is then derived from the difference between the outputs of the 
photodetectors. 



Using TTt(X) = <Jt(X) I <J t {I) and the Ito rule, we find the 
following homodyne quantum filter 

dTrt(X) = n t (£(X))dt + 

a (lit ( sin( K F, )X + X sm( K F z )) - 2n t ( sin{nF z )) ir t (X)) 
x (dY t - 2an t ( shi(kF z )) dt) . 

(20) 

The process dY t — 2onr t ( sm(nF z )^jdt is again called the in- 
novations or the innovating martingale. It can be shown 
[To! Ihij that the innovations are a continuous martingale 
with respect to the filtration y t and the measure induced 
by P. It follows from Levy's theorem that the innovations 
for the homodyne detection setup form a Wiener process. 

VI. STRONG DRIVING, WEAK COUPLING 

Define the measurement strength as the product M = 
a 2 K 2 . In a typical experimental setting a will be very 
large (strong driving) and k will be very small (weak 
coupling). The idea in this section will be to exaggerate 
this by taking the limit a — > oo while keeping the product 
AI = a 2 n 2 constant. 

Let us introduce the following scaled sum and differ- 
ence processes 



aa t (ain{KF z )X + Xsm(K,F z ))dY t . 



(19) 



Yt + Y? 



y; 



Yr 



(21) 
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Note that we scaled the sum by a 2 and the difference 
by a. We will see that with these scalings we get finite 
output processes in the limit. In practice the scalings 
are determined by the experiment, i.e. they are chosen in 
such a way that the photocurrents nicely fill the scales on 
the read out devices. We are interested in the statistics 
of the processes Y t + and Y t ~ . Therefore, following [24| . 
we introduce their characteristic functionals 



$ + (fc,t) = P ( cxp \-i J k{s)dY s 

rr* ( ■ f* fc ( fi ) 
^ CXP Wo ^ 

$~(M) = F ( CX P (-i J Hs)dY s 



(dAr + dAf))[/ t ), 



(dAf + dAf ) )U t ), 



= P [U; exp -i 



where k is an arbitrary function in L 2 (M. + ). The char- 
acteristic functionals <I> + and <I>~ faithfully encode the 
complete statistics of the processes Y t + and Y t ~ . 

Using the quantum Ito rule and the fact that vacuum 
expectations of stochastic integrals are zero, we find the 
following differential equation for <&+(£;,£) 



~~dt 



(k,t) 



= a 2 | exp [ -i-Q- | - 1 



*+(M), 



In the limit a to infinity (while keeping M = a 2 K 2 con- 
stant), we therefore obtain 



$+(fc,t) = exp (— i J k(s)ds 



This is the characteristic functional of the deterministic 
time process t. In short, as a tends to infinity, dY^~ tends 
to dt. 

To calculate $~(k,t), define for all atomic gas opera- 
tors X 



$-(X,k,t) 



UfX cxp ( —i 



k(s) 



(dA x s y + dAf ) U t 



Note that s) = $"(/, k, s). Using the quantum Ito 

rule and the fact that vacuum expectations of stochastic 
integrals are zero, we find the following system of differ- 
ential equations (n > 0) 



(a n sm n (2KF z ),k, 
a 2 (cos - l)^(a n sm n (2 K F z ),k, 

iasin $~(a n+1 sin" +1 (2 K F z ),fc,t 



Note that although the atomic gas system might be very 
high dimensional, the dimension is finite. That means 



that the above system of differential equation is closed 
and consists only of a finite number of equations. In the 
limit a to infinity (while keeping M = a 2 K 2 constant), we 
obtain the following finite system of coupled differential 
equations (n > 0) 



d<£- 
~dt 



(VMF z ) n ,k,t) =—U-$-((VMF z ) n ,k,t) 
2ik(t)<f>~ ((yfMF z ) n+1 ,k,t 



(22) 



In principle we could now try to solve this system of equa- 
tions. However, instead of finding an explicit solution, 
let us compare this with the statistics of the homodync 
observations Y t defined in Eq. (TTT)) . In analogy to the 
discussion above, we define for all atomic gas operators 
X 

${X,k,t) = 

Using the quantum Ito rule and the fact that vacuum 
expectations of stochastic integrals are zero, we find the 
following system of differential equations (n > 0) 



— (a" sin" (KF z ),k,t 



k{tf 



$ (a n sin n (KF z ),k,t) - 



2ik{t)<p(a n+1 sin" +1 (nF z ),k, 



Taking the limit a — ► oo (while M = a 2 K 2 is held con- 
stant) then again leads to the system of differential equa- 
tions Eq. (|22j) . Therefore we conclude that in the limit 
the processes Y t ~ and Y t have exactly the same statistics! 
This means that from the point of view of statistical in- 
ference of the atomic gas system from the observations, 
the balanced polarimctry experiment and the y-channcl 
homodyne detection experiment are equivalent. 

Rearranging terms, we can write the linear quantum 
filtering equation Eq. (fT5|) as 



da t (X) = a 2 <j t (sm(KF z )X sin(KF z )) 



a t (cos(nF z )X cos(kF z )) - o t (X) \dY t + 



<T t (cos(KF z )X sin(KF z )) + 



a t (cos(nF z )X sin{ K F z )) \dY{ 



Writing Y t for the limit process of Y t , and taking the 
limit of the above equation, we obtain the following linear 
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quantum filtering equation 

da t (X) = a t (Z(X))dt + VMa t (F z X + XF z )dY t , (23) 
where 

C(X) = M(F,XF, - \{F 2 X + XF 2 )). 

Moreover, we obtain the following normalized quantum 
filter 

dw t (X) = TT t (C(X))dt + VAl(n t (F z X + XF Z ) - 

(24) 

27r t (F 2 )7r t (X)) (dY t - 2VMv t {F z )dt) . 

Since dY t — 2\/M7r t (F z )dt is a continuous martingale 
it follows from Levy's theorem that it is a 
Wiener process. That is, we find that dYt = dWt + 
2\[M'Kt(F z )dt, with W t a Wiener process. 

Furthermore, note that if we start from Eq. (fT!?|) . tak- 
ing a to infinity while M = a 2 K 2 is held constant, then 
we also obtain the linear filter Eq. (|2"3"|) . Likewise, the 
homodyne filter Eq. ([2H| converges to the filter in Eq. 
([2"4")l when a is taken to infinity while M = a 2 K 2 is held 
constant. 



VII. DECOUPLING THE X-CHANNEL 



Now, if we replace k by \[M ja and take the limit a to 
infinity, then the time evolution satisfies the following 
QSDE 

dU t = | VMF^e^dAH - e^dAf) - ^-F 2 dt^U t . 

That is, the x-channel has been decoupled from the in- 
teraction, see also [1, HU . 

In a similar way we easily see that in the strong driving, 
weak coupling limit we have dY^ = dt and for Y t , the 
limit of Y t ~ , we obtain 

1 1 — 

a 

_ U* t {ae~ 1 ^ A\ + ae^ApUt (25) 
a 

= U*(e- i<t,t A y t + e^ADUt. 

This shows once more the equivalence of the balanced 
polarimetry experiment and the y-channel homodyne de- 
tection experiment. It is easy to see that the charac- 
teristic functional of the process Y satisfies the set of 
coupled differential equations of Eq. (|2"2"1) . Moreover, 
after decoupling the x-channel the system is given by 
jt(X) = U* t XU t and the observations by Eq. ((25j). Fol- 
lowing essentially the same steps as in Section |V1 we 
again obtain Eqs. (j23|) and (|24|) as the linear and nor- 
malized filters, respectively. 



Let us give a brief formal discussion to show what hap- 
pens in the strong driving, weak coupling limit. As a 
increases and k = \lMja decreases, the relative effect 
of the atoms on the x-polarizcd channel also decreases. 
Therefore, we can reasonably expect that the x-channcl 
remains in a coherent state. Instead of working with re- 
spect to the state P = p®<j> we now work with respect 
to the state 

q = p® (</>*(/), ■ ru)), 

with f{t) = ae*. Note that this means that the y- 
channcl is still in the vacuum state. Working with respect 
to the coherent state on the x-channel means that the 
time evolution is given by Eq. (JJ) . Formally we can write 
for /3, 7 e {x,y} 

dAf 7 = a f l*ajdt, dA p t = afdt, dAf* = af*dt, 

and since for large a and small k the x-channcl is approx- 
imately in the coherent state ip{f)> we can replace by 
ae"^' and af* by ae -2 ^' . This means that we obtain for 
large a and small k 

dU? = { ( cos(kF z ) - 1) (a 2 dt + dAf ) - 

sin( K F z )a(e-^dA y t - e^dAf )}lJ°. 



VIII. DISCUSSION 

We have provided a quantum stochastic model Eq. © 
to describe recent polarimetry experiments in which po- 
larized laser light interacts with an atomic gas via the 
Faraday interaction. In our description the gauge pro- 
cess plays a prominent role. It represents the scattering 
between different channels in the field and it provides us 
with counting processes that can be observed. As in [2l| . 
our quantum stochastic model presents a novel applica- 
tion to quantum optics of the gauge terms in a QSDE. 

Once we set the model, we derived quantum filter- 
ing equations for balanced polarimetry and homodyne 
detection experiments, studied the statistics of output 
processes and obtained filters in the strong driving/weak 
coupling limit. Our results in the limit confirm the ad 
hoc filter for the balanced polarimetry experiment that 
has already been in use in the literature [IJ. Moreover, we 
showed that from the point of view of statistical inference 
the balanced polarimetry experiment and the homodyne 
detection experiment are equivalent. 

Using formal arguments we have seen that in the strong 
driving, weak coupling limit the a;-channel decouples 
from the description. Rigorous results on this decoupling 
are still to be obtained. 

Having an underlying model from which rigorous 
derivations can depart, is likely to be advantageous in 
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future further investigations. In particular, combi ning 
the model presented in this paper with the results in [261 ] 
could prove useful for investigating the situation where 
the laser beam passes through the gas multiple times 
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